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Abstract. We construct solutions, for small values of G and angular fre- 
quency f2, of special relativistic scalar gravity coupled to ideally elastic mat- 
ter which have helical but no stationary or axial symmetry. They correspond 
to a body without any symmetries in steady rotation around one of its axes 
of inertia, or two bodies moving on a circle around their center of gravity. 
Our construction is rigorous, but modulo an unproved conjecture on the 
differentiability of a certain functional. 
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1. Introduction 

In Newton's theory of gravity there are solutions describing two bodies which 
move at constant angular velocity on a circle around their common center of 
gravity. The existence of such solutions for fluid bodies was demonstrated by 
Lichtenstein 1922 [LI]. In Einstein's theory we do not expect such solutions to 
exist, because such a system will emit gravitational radiation, and so the two 
bodies will spiral towards each other. 

However in 1992 Detweiler and Blackburn |DEj conjectured that such so- 
lutions should exist also in general relativity, provided there is just the right 
amount of incoming gravitational radiation to keep the bodies on their circle. J. 
Friedman, R.Price and coworkers ( |YOj . |BEj and references therein), followed 
this idea and analysed various model systems mathematically and numerically. 
There are many conceptual and technical problems involved. From the Schild 
solution |SCHj of electromagnetism, which describes two point charges moving 
on a circle, it becomes clear that such a system has infinite energy because the 
required amount of incoming (and outgoing) radiation is unbounded. So, in GR 
we expect that the ADM mass will be infinite. Consequently the system will 
not by asymptotically flat in the sense of any of the known definitions, which 
in turn makes it difficult to even define such a system. We want the spacetime 
to admit a Killing vector which behaves like dt + ^d^f,. The meaning of this, 
however, is unclear if there is no asymptotic symmetry group. There are ways 
out, Friedman proposes a definition in |Y0] . 

There is another subtlety: in Minkowski space the helical Killing vector dt + i^d^ 
is timelike near the center and spacelike near infinity. Hence, the metric on the 
quotient of the Killing vector changes signature. It is completely unclear, how 
to deal with the Einstein equations in such a situation. 



Date: May 18, 2009 File:Hehcal.rot.body7.tex . 

* Supported in part by Fonds zur Forderung der Wissenschaftlichen Forschung project no. 
P20414-N16. 

1 



2 



R. BEIG AND B. SCHMIDT 



All numerical model systems treated so far do not treat the 'bodies' dynami- 
cally. In this paper we consider a special relativistic theory in which the bodies 
are composed of some elastic material and gravity is modeled by a scalar field. 
The main point of the paper is to demonstrate that solutions do exist in which 
the incoming radiation keeps the bodies on their circle. In special relativity 
there is no problem to define such systems. We just want the solution to be in- 
variant under the Killing vector dt + ^d^f,. The scalar gravitational field satisfies 
the wave equation in Minkowski space. 

There is a case which is technically simpler than the 2-body system. This is 
a rigidly rotating 'tri-axial body'. By this we mean that the body does not have 
to have any symmetries. Again, in Newtons theory such solutions exist |BS3j . 
In GR, however, the body should spin down because it radiates gravitational 
waves if it is not axisymmetric with respect to the axis of the rotation. 

The plan of the paper is as follows. In Section 2 we describe the theory we 
are using: special relativistic elasticity formulated as a Lagrangian field theory 
|BS1| minimally coupled to a version of scalar gravity. Section 3 discusses 
helical symmetry and derives the equations for which we want to show existence. 
In Section 4 we consider a triaxial body in helical motion and describe 'its 
gravitational field' according to scalar gravity. We try to find a solution with 
the right amount of incoming radiation by taking as gravitational field (1/2 
of) the sum of the retarded and advanced solution of a given source motion. 
With any other combination our method would break down. We want to show 
existence using the implicit function theorem. Because the linearized operator 
has a nontrivial range we first solve a 'projected system' first and only then the 
full system. 

In Section 5 we discuss in outline the helical 2-body problem for which exis- 
tence can be established by essentially the same method. We were not able to 
show that no solution exists if we take just the retarded solution of the wave 
equation. 

There is a gap in our existence proof, because we do not demonstrate that the 
'gravitational force' depends differentiably on the configuration, in the function 
spaces we are using for elasticity. We consider this as a technicality. If this point 
can be cleaned up, we are confident that even for complicated non linear systems 
such as GR the idea that incoming radiation can balance outgoing radiation 
works. Finally, we use elastic bodies only because in this case we understand 
the boundary value problem. A corresponding result should certainly also be 
true for fluids. 



2. Scalar gravity coupled to elasticity 

Consider a Lagrangian field theory on Minkowski space M space with metric 
^ij,u = diag(— 1, 1, 1, 1) with dynamical variables given by a scalar field V{x^^) 
descibing gravity and and the fields f^{x^) describing elasticity in Special Rel- 
ativity [BSlj . which are viewed as maps M ^ B with B, the so-called material 
space or body, being a bounded domain in with smooth boundary dB. The 
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action is 

S = IJ V^''V^.VuV^ d^x + 47rG j pF{V)y/^ d^x. (2.1) 

We assume that B is endowed with an Euchdean metric 5ab^ where A,B = 
1,2,3 and that the scalar p defining the elastic material depends only on the 
principal invariants of the matrix H^b = 5bc^ where H^^ = f^^^f^^uV^'^- 
This means, in standard language, that we consider isotropic materials. The 
energy momentum tensor of the system is: 

= yi^^,u - \ri^.uV''\aV,x + ^ttGT^^F, (2.2 
with r^jy the elasticity energy momentum tensor from |BS1] . 

Tfj.u = pu^Uy - a^y , (2.3 
where is the future-pointing vector field given by 

/^,^n^ = 0, r^^.u^^u" = -I (2.4 

The existence and uniqueness of is a regularity condition on the fields 
we are considering. The quantity a^^ is given by 

— '^^ ^9 ^ 

where 

p = ne , (2.6 

and n is defined by 

eABcf^Ax)f''Ax)fA^) = «(^) (2-7 
The equation for V is 

UV - AirGpF' = 0. (2.8 
To obtain a linear equation we choose 

F{V) = l + V (2.9 

We have 

V^r*°* = 0. (2.10 

The equations (|2.2I2.1U|) imply 

= pV^, + T^^V^ + (1 + V){pu^,u^ - cj/),^ (2.11 

Or {h^- = b^^ + u^u-) 

= \p h^'' - + (1 + y) {pu^v" - a^")^, (2.12 

Note that in the calculation of a^/ ,p in terms of the the first term in p does 
not contribute. Furthermore, both the 'force term' (involving Vy) in (j2.12p and 
the remaining term is orthogonal to (for the latter this is e.g. shown in 

mm- 

We have a PDE system for V, f"^. This has a Newtonian limit which is treated 
in Sect.5 of [BS3] . 
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3. Helical motion 

We assume that the material flow is parallel to the helical Killing vector 
dt + ^d^, i.e. that 

f^^,,{dt + nd^Y = (3.1) 

and that V is invariant under dt + ^29^, i.e. 

v^{dt + nd^r = o (3.2) 

For concreteness we choose coordinates (r, y) , so that r = coincides with the 
hyperplane t = with y Euclidean coordinates thereon and dt + fic?^ = dr- 
Explicitly the transformation is given by 

= y^ cos Qt — y"^ sin fir, = y^ sin fir + cos = y^, t = t (3.3) 

In these coordinates the configuration can be written as functions /"^ of y, 
and similarly for V. By slight abuse of notation we will omit the bar in what 
follows. Now Eq. (j3.ip implies that the quantities n, p and a^'^ are all also invari- 
ant under dt + 0,8^ and (T^^{dt + ^d^Y — 0- Using the orthogonality pointed 
out at the end of the previous section and the expression for the projection of 
(TfY ,v on the quotient calculated in the paper |BS3j . we obtain 

= [phj' - cj/] DiV + {l + V) [-e-^Di{e^a/) + pDjU] , (3.4) 

where — e^^ is the norm of the helical Killing vector given by 

e2^ = l-0V, = (x^)^ + {x^)\ r^n^<l (3.5) 

hij is the metric on the quotient of dt + $19^ and Di the covariant derivative 
with respect to htj. Note that, for $7 = 0, hij equals 6ij. Note also that now 
jjAB ^ jA .fB .^ij_ rewrite ([33]) as 

A [(1 + V) e^a/] -pe^ [DjV + (1 + F) DjU] = (3.6) 

The second term and third term in ()3.6p is respectively the gravitational and 
centrifugal force. We write ()3.6p as 

= A [(1 + V) e^a/] + Wj + Zj , (3.7) 

and 

-Wj = e^pDjV (3.8) 

with 

- Zj = e^p{l + V)DjU , (3.9) 

and 

DjU = n'^ {l-n'^r^y^' ^Djr^ (3.10) 
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4. Triaxial Rotating body 

The first case we want to consider is a triaxial rotating body coupled to 
scalar gravity. We repeat that 'tri-axial' means that the set B need not have 
any symmetries. Let us denote the corresponding relaxed configuration by $. 
We take for the gravitational field the symmetric (i.e. one half the sum of 
the retarded and advanced) solution of the scalar wave equation. Note that 
this makes sense globally in Minkowski space. But we will only make use of 
that solution inside the support of the body, where r^O^ < 1. For a source p 
invariant under helical motion we have from the Appendix that 

V{y) = G J H{y,y')p{y')y/-e{y')Hy')dy', (4.1) 

where H(y,y') is given by (16.12p . Clearly this function is even in Q. Expanding 
V m Q only even powers occur and we have 

V = GV = V^ + + • • • = G{V^ + Q'^V^ + ...) (4.2) 

where is the Poisson integral of the source. 

4.1. Solution of the projected system. The fundamental equation (with V 
still chosen half the sum of retarded and advanced potential defined by p) is 

= Di [(1 + V) e^a/] + GWj + n^Zj (4.3) 

with 

Wj = GWj (4.4) 

and 

Zj = n'^Wj (4.5) 

Equation (|4.3p forms a quasilinear second-order system of partial differential 
equations for the functions f"^, which has to be solved subject to the (Neumann- 
type) boundary condition 

= (4.6) 

Due to the presence of a free boundary we will also need the material form of 
(j4.3p . where the dependent variable f{y) is replaced by its inverse ^{X) with 
X eB. Then, with Wj = n-^Wj{(l){X))) and Zj = n-^Zj{(l){X))), there results 

= Fj{<^, n^, G) = Va [(1 + GV) e^a^j] + GWj + ft'^Zj (4.7) 

where a^j is the 'first Piola stress' given by a^j = n^^ f^ia^ j. The boundary 
condition now takes the form 

(y^inA\dB = (4.8) 

We have two parameters in the problem, namely and G. 
Now one has to define an appropriate space for ^{X), which will be a neigh- 
bourhood of the identity map in a suitable Sobolev space, we refer for details 
to the paper |ABSj . Let $ be a relaxed, i.e. stressfree configuration, which 
we take to be the identity map. We also take B to be contained in the set 
[(X^)^ -I- {X'^)'^]Q? < 1, whence, for deformations ^{X) sufficiently close to the 
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identity, the set ^{B) is inside the hght cyhnder [(y^)^ + < 1. 

We next require that the function e{H^^) in ()2.6p satisfies 



e(5^^) = e = const > 



= (4.9) 



qhAB 

Then we have a solution Fj{^^ 0, 0) =0 of (|4.7p . Assume that furthermore 

= A 5ab5cd + 2/i 6c[A^B)D (4.10) 



( \ 

V dH^BQHCD 



HEF= SEF 

with the constants fi, A obeying i^l > 0, 3A + 2/_f > 0. The linearization of Fj at 
($, 0, 0) is the standard operator of flat-space linear elasticity on B, which has 
to be considered together with the linearized form of the boundary condition 
()4.8p . It is well-known (see e.g. |ABS] ) that this is an elliptic operator with 
finite dimensional kernel and range. The latter is given by fields li{X), such 
that Jjg i^^iX) 5^ A k{X) d?X = for all Euclidean Killing vectors r/"^ on B. One 
can like in [ABSj define a projection P onto the range of the linearized operator 
and study the projected system 

= PFi($, G) = P \Va [(1 + GV) e^cT^i] + GWi + O^Z^} (4.11) 

Assuming differentiability of this nonlinear map we can apply the implicit func- 
tion theorem. Differentiability is standard for the Zi - term (see |BS3] ) . For the 
Wi - term, which contains the retarded -|- advanced scalar field, we have to leave 
this as a conjecture. The corresponding result for the Newtonian gravitational 
field has been proved in [BS2j . 

The kernel of the linearized operator consists of Killing vectors of flat Euclidean 
space. We restrict the deformations <^ to ones of the form 

¥{X^, G, n^,a\ V^) = b^AX^ + a* + }^^hAX^ + l> (4.12) 

where a* and V"^ are constants with = b^^^^ and 4> is in a complement of the 
kernel which contains ^. Thus we obtain a unique solution of the projected 
equations for each value of a* , b^^ defining an infinitesimal motion. Hence we 
know ^{X'^,G,Q'^,a\b'^). 

4.2. Equilibration. To obtain a solution of the full equations we will construct 
a family {a'{G, Q"^), b'''{G, n"^)) such that 

^\X^, G, n^, {a\G, 1^2), 0^))) (4.13) 

satisfies the equation and the boundary condition. We will determine this family 
by solving (I — ¥)Fi = 0. This condition is equivalent to 

= %)(G,172,C) = / ±. {Va [(1 + GV) e^a^i] + GWi + ^^Z,} (fX 
Jb 

(4.14) 

for all Killing vectors ^^^^ with a = 1, ...6 of the flat metric on the body, 

C = {a\b''') and $*(X^, G, Q^, G) is inserted in the integrand. Note that 
^1^^ = a'* -|- b'^j-y^ has to be composed with <I>, so the rotational terms depend 

on ^^{X"^, G, O^, G). To have only one parameter we assume u;^ = kG with k 
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a positive constant and, with this in mind, define the "normahzed force map" 
as 

A^(,)(G, C) = G-^ [ 4) {Va [(1 + GV) e^a^i] + GW, + V?Z,] S'X (4.15) 

JB 

and replace (|4.14p by 

iV(„)(G,C)=0 (4.16) 

This is well defined because both the forces and of have a factor G. 
In analyzing the 'equilibration condition' (|4.16|) we show, first of all, that (|4.16p 
is indeed satisfied for arbitrary G, when ^(^^ is dyi or ip'dyX — y^dyi or a constant 
linear combination thereof. To prove this, it is simplest to look at the spatial 
version of ()4.16p . namely (see ()4.3p ) 

G~^ I CL) I A [(1 + V) e^a/l + GWj + n^zA Vhcfy = (4.17) 
J'S'iB) ^ ' ^ ^ 

Now, the first term in (|4.17p is zero: use integration by parts, the Killing 
equation for each of these two vectors and the vanishing of the boundary term 
by virtue of (|4.6p . For the third contribution already the integrand is zero, 
again by the symmetries of U. Finally, the second term is zero: this amounts to 
the statement, proven in the Appendix, that, for a body which rotates rigidly 
around the - axis, the - component of both the force and torque due to 
its own gravitational field is zero, provided this field is given by the symmetric 
solution of the wave equation (12. 8p . 
We thus only consider 

iV(„,)(G,C') = (4.18) 
where C' = (o^, o^, a"^ = 0, 6^^ = 0, 6^'^, 6^^) and ^(^f') to be the collection 

U(i),^(2),C(i3)>^{23)} = {d^i,dy2,y^dyi -y^dy3,y^dy2 -y^dyi}. (4.19) 

Suppose A'"(q/)(0,0) = and -^-^(0,0) is invertible. Then, by the (finite- 
dimensional) implicit function theorem, for sufficiently small G there exists a 
function G'{G) such that N(^^i){G,G' {G)) = 0, and we have a solution of our 
problem. 

Calculating iV(0,0) is the same as The first term in the integrand gives 
dA^crf which is equilibrated as was discussed in Eq.(4.2) of |ABSj . So we have 
only to consider the force terms. 

N^^^,) (G, C') = ( {Wi + kZ,) d?X . (4.20) 

We obtain 

iV{a')(0'O=- / 
JB 

(4.21) 

where ^ = b\X^ + a* + V-^bkA^^ with = 6^^ = q. The - term is 
simply the Newtonian gravitational self- force (resp. self-torque), and so its 
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contribution to the integral vanishes (see e.g. example (i) in the Appendix). 
We are thus left with 



where 



iV(i)(0,C") 
iV{2)(0,C') 
%3)(0,C") = 

iV(23)(0,C") = 



(FX 



- [ {X^+a^ + b^^X^) d^X 
JB 

JB 

j X^ (X^ + +b^^X^)d^X 
JB 

[ X^ {X^ + a'^ + b'^^X^)d^X 
JB 



(4.22) 

(4.23) 
(4.24) 
(4.25) 
(4.26) 



From here we infer that A''(-„/)(0, 0) is zero iff (i) the rotation axis in the reference 
estate goes through the center of mass and (ii) coincides with one of the axes of 
inertia. In fact, we assume that the center of mass is at the origin, rather than 



just (i). With these assumptions minus the matrix ' (0, 0) is given by 



/ V 

V 



\ 



where 



V 



d^X 






W 


w 



\ 




w / 



(4.27) 



dN, 



So ' (0, 0) is invertible, and our argument is complete. 

5. Two BODIES IN CIRCULAR MOTION 

In this section we consider two identical bodies in circular motion around 
their common center. To simplify matters we assume that the support of the 
matter in the relaxed configuration consists of two spherical balls with centers 
at {—L, 0, 0) and (L, 0, 0) on the quotient metric. We have as discrete symmetry 
the reflections {y^,y'^,y^) {y^,-y'^,y^) and {y^,y'^,y^) {y^ ,y'^ , -y^). For 
an isotropic stored energy functions we know that the solution will also have 
this symmetry. Therefore we consider only configurations with the property 

^>l(X^X^X3) = ^^{X^,-X'^,X^) = $^(X\X2,-X3) (5.1) 



^>l(x^x^x3) = -<^^{x^ 



-x^,x'-^ 



^>^(x^x^-x' 



^>^(x^x^x2) = ^^{x^,-x'^,x^) = -$^(x^x^-x3) 

One is now able to show the following: 

1.) For such configurations the gravitational field "ad+ret" inherits this 
symmetry. So does the centrifugal force. 



2.) For an isotropic stored energy also the stress tensor and its divergence 
inherit this symmetry because it is an isometry of the quotient metric. 
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3. ) All these properties imply that all equilibration integrals, with the ex- 
ception of that for the - translation Killing vector, are satisfied for any 
configuration with these symmetries. 

4. ) If we assume further that y^ — > —y^, which exchanges the two bodies, is 
a symmetry we have only to consider the elastic equations for one body. 

5. ) As usual, we first solve the projected equations with a parameter a^, 
the translation in the y^-direction, and have then to study the bifurcation 
equation N{G,a^) = 0. It turns out that N{0,a^) = is the same as in the 
Newtonian problem which we solved in |BS4j . We find a reference configuration 
with A^(0, 0) = and ^(0, 0) = and have thus solved the 2- body problem. 

6. Appendix 

In this Appendix we prove a result which contains as a special case the state- 
ment that the gravitational contribution to the equilibration integrals (|4.15p 
vanishes, when the Killing vector is taken to be d^s or (9^ = —x^d^i + x^d^2. 
Let p{t, x) be a source and 

V(..)=/g(*-*';.-.')p(*',.V*VV (0.1) 

be the field, where G is 'some unique' Green function (in fact: distribution) of 
□<1> = p. By 'some unique Green function' we mean a Green function sharing 
the symmetry of the background (i.e. Minkowski space in our case), such as the 
retarded or advanced Green function. We assume that p has compact support 
in space. 

Let us next suppose that the source p is invariant under the flow of some 
timelike Killing vector ^^d^. Using coordinates comoving with ^, i.e. (t; y) so 
that i^d^jL = dr, Eq.(l) takes the form 

V{t, y) = j Hiy, y')p[y'U -e{y')h{y') , (6.2) 

where H{y,y') = f G{t — t' ]y,y')dT' , .^^ = g^j.ui^i'^ and h is the determinant 
of the metric on the quotient of (The square root in (2) is \/—g in the 
comoving system.) Clearly V is actually independent of r. 
Let now be a spacelike Killing field commuting with ^, i.e. projecting to 
^ = rf{y)di on the quotient. Consider the expression given by 

Fn = j P{y) Wdi H{y, y')] p{y') ^-^^{y)h{y) V-C'(y')%') d^v' (6-3) 

This is up to sign the equilibration integral Jj-i^jg^ir]''WiVhd^y. Since the 
(unique) Green function shares the symmetry of the background there has to 
hold (operate on both arguments with the symmetry and linearize) 



W{y)di + /iy')d,]Hiy,y') = 



(6.4) 
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(Note this relation holds irrespective of whether we take the retarded, advanced 
or any other combination.) Furthermore rfdi is a Killing vector in the quotient 
space, so that in particular 

diiVhri') = (6.5) 

Also there holds 

^'di f = (6.6) 
We use (j6. 416.516. 6p and integration by parts in y' to obtain (the two minus-signs 
arising in the process cancel) 

F, = J p{y) H{y,y') W\y')d, p[y')W-e{y) Ky) ^-^iv') Ky') d^d^y' 

(6.7) 

But we can also perform partial integration in the y-variable to find that 



Pv = -J [v\y)^^ p{y)] Hiy, y') p{y')^-e{y) Kv) V-e{y') Ky') d'y A' 

(6.8) 

We now assume that the Green function is the symmetric one, so that G is 
invariant under the joint interchange of {t,t') and [y,y'). It follows that 
H{y,y') = H{y',y). Whence, by (I6.3|6.8|) . we find that F,, is zero. We il- 
lustrate the previous result by two examples: 

(i) static symmetry: Here ^ is given by ^^dfj_ = dt, and the, say retarded. Green 
function in adapted coordinates takes the well-known form {t = t, y = x): 

,^G..,^r-r';y,y') = '^^-;'-\-''\^ (6.9) 

\y-y\ 

Furthermore 

A'KH{y,y') = —^ (6.10) 

\y - y \ 

i.e. H{y, y') is the Poisson kernel. (This result is of course independent on 
which combination of the retarded and advanced Green function has been cho- 
sen.) The quotient metric h is the Euclidean one and spacetime Killing vectors 
projecting to the quotient in this case are simply all Euclidean ones in y - space. 
We have thus recovered the well-known Newtonian result, that the force and 
torque on a static body due to its own gravitational field are zero. 

(ii) Helical symmetry (see jBHSj ): Here ^ is given by ^^d^ = dt + ^d^, and the 
retarded Green function takes the form (t = t, p = cj) — Qt, z = x^) 

47rGrct(T - T';y,y') = 



6 (t -t' - Y^r2 + r'2 - 2rr' cos[p - p' + ^{t - r')] + {z - z'Y 



(6.11) 



^j,2 _|_ _ 2rr' cos[/.i — p' + Q.(t — t')] + (z — z')'^ 
whence (assuming r'^Q,^ < 1) 

A7rH{y,y') = U— J-^ --— J-^ -), (6.12) 

2 \a^{p — p,r,r , z — z' ) a-[p — p,r,r, z — z' ) J 

with cr±{p, r, r' , z) being implicitly given by 

a± = ± y^r2 + r'2 - 2rr' cos{p + na±) + z^ (6.13) 
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(If we had not chosen the symmetric Green function the quantity correspond- 
ing to H(y,y') would have been different, in particular not symmetric.) The 
quotient metric h is given by 

hij dx'dx^ = dr^ H d^^ + dz^ (6.14) 

1 — ii^r^ 

and spacetime Killing vectors projecting to the quotient are given by d^z = dz 
and (?0 = dfj_. Applying the main statement of this Appendix to these two 
Killing vectors yields the result concerning the second term in Eq. (j4.17p . 

Acknowledgement. Jiirgen Ehlers has greatly influenced our scientific think- 
ing. We dedicate this work to his memory. 
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